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Abstract 

The nucleation and/or growth of cracks in elastic-brittle solids has been recently 
described in [14] in terms of a special class of measures and with a variational 
technique requiring the minimization of a certain energy over classes of bodies. 
Here, the physical foundations of the theory and the basic ideas leading to it are 
described and commented further on. A view on certain possible developments and 
shifts toward different settings is also given. This article has expository character. 
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When the adjective 'variational' is attributed to models of some physical phe- 
nomena, impulsively one may think that the physical situation under scrutiny 
is conservative. Variational means, in fact, that one is managing some func- 
tional, essentially an energy, and is asking that it attains its minimum on some 
function space. The situation is well known and is typical, for example, of the 
determination of ground states of elastic bodies under prescribed boundary 
conditions. The energy is defined on the reference place B occupied by the 
body in the ambient space, a place fixed once and for all. Its minimum is 
required to be attained over a class of one-to-one orientation preserving maps. 

Variational descriptions of the nucleation and the evolution of defects like 
cracks however exist. At least in principle, the question of their physical ap- 
propriateness can be posed, due to the dissipative character of the phenomena 
they are referred to. Of course, one can say that the quest of something to 
be minimized has rather instrumental character because in nature some eco- 
nomic principle always appears somewhere and the calculus of variations is 
at a stage of development that it can be desirable to be under conditions of 
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using it. This point of view can be an answer. However, the answer can be 
even deeper. 

When one think of the nucleation of defects, in fact, one manages a mutant 
body. Mutations occur, in fact, in the gross shape of the body at the contin- 
uum level, and they can be naturally pictured as mutations of the reference 
place B which is now not fixed once and for all as for example in the standard 
formulation of elasticity. In principle, one can imagine of having at disposal a 
class of possible bodies occupying B, every member of the class differing from 
the others by the defect pattern. Once boundary conditions are prescribed, 
a variational principle can select the resulting body within a class. This way, 
a variational approach to the defect nucleation - the latter imposed by the 
interplay between boundary condition and nature of the material - is physi- 
cally significant when it involves a minimization process of the free energy (for 
example) over an entire class of possible bodies. This point of view is the one 
adopted here and is also the answer to the question of physical appropriate- 
ness of some variational approaches. The energy dissipated in the nucleation 
(or growth) process is individuated in the gap arising from the body in the 
original shape - the one existing before the assignment of boundary conditions 
- and the actual shape (i.e. the body plus the defect pattern) obtained by the 
minimization procedure. Of course the explicit evaluation of this gap can be 
arduous. 

This point of view interprets the physical aspects of two existing variational 
models of fracture mechanics. Both models are recalled below. The attention is 
focused on the second one in which cracks are described through appropriate 
measures, the so-called curvature varifolds. The physical significance of the 
approach is discussed here. It is shown also that the point of view can be 
adopted in other circumstances dealing with the nucleation of defects stratified 
over manifolds with different dimensions. 
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Under the suggestion of the pioneering Griffith's approach to fracture in 
brittle-elastic solids, a variational model in fracture mechanics has been pro- 
posed in [10J. It is based on a requirement of minimality of the overall energy 
(£ which accounts for the macroscopic deformation and the possible presence 
of cracks, and is defined by 

£(C,u) := / e(x,Du(x)) dx + [ (p dV. 2 , 

J B JC 
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where B is the reg region occupied in the three-dimensional ambient 

space by the body, Du (x) the spatial derivative evaluated at x G B of the 
differentiable transplacement map x i — > u (x) G M. 3 - the map defining the 
actual (deformed) place u (B) - C the representation in B of a surface-like 
crack occurring in u (B) , e the elastic energy, <f> a constant surface energy, ctH 2 
the two-dimensional Hausdorff measure. The map x i — )• u is also assumed, as 
usual, to be one-to-one outside C, orientation preserving (i.e. det Du (x) > 0), 
and such that the global invertibility condition 

/ / (x, u (x)) det Du (x) dx < sup / (x, z) dz, 

holds for all / G C^°(B x IR 3 ), a condition allowing frictionless self-contact 
of the boundary while still preventing self-penetration (for details on this 
last condition see [17]). In other words, x i — > u is an orientation preserving 
homeomorphism outside the subset of C containing its jump set. 

The requirement is then that at each instant t G [0, t\ of a cracking process 
the pair (C, u) (t) realize a minimum of the global energy (£ with C an admis- 
sible crack. Admissibility is intended in the sense that C is a rectifiable set 
(specifically the image of a countable number of Lipschitz maps) with zero 
volume measure. The interval of time is then discretized and minimality is 
required at time steps. Minimizers are pairs (C,u): the minimum problem has 
two variables. 

Notwithstanding simplicity and elegance of this model, the evaluation of min- 
ima of the energy at each time step involves a number of analytical problems. 
The essential difficulty arises in controlling in three dimensions minimizing se- 
quences of surfaces leading to the possible actual crack, or better to its picture 
C in the reference place. However, when the entire crack is open, C coincides 
with the jump set of the transplacement field x i — > u(x). The convenient 
simplification of identifying cracks with the jump sets of displacement fields 
has been then adopted. Bounded variation (BV) or special bounded variation 
(SBV) functions have been then involved as candidates to be minimizers of 
the energy considered as a functional of the sole u. The energy is then con- 
sidered as the one of a simple Cauchy's body^l and minimizers are sought 
in a space of maps including candidates to be reasonable descriptors of the 

1 The type of regularity is specified later. 

2 In continuum mechanics a body is considered in primitive sense as a set with 
elements called the material elements and identified even vaguely with molecular 
or atomic aggregates. The selection of material elements, also called representative 
volume elements, is matter of modelling. Once even a rough idea of them is formu- 
lated - in a sense one is specifying what are the peculiar physical aspects of the 
material texture - the essential step is to furnish geometrical structure to the body 
which would be otherwise just an abstract set. The representation of the interactions 
is then straighforward, dual in the sense of power. The geometrical representation 
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elastic-brittle behavior. Remind that the choice of function spaces where one 
researches minimizers has constitutive nature. A difficulty has been however 
encountered: theorems allowing the selection of fields with discontinuity sets 
describing reasonable (physically significant) crack patterns seem to be not 
available at least in the current literature (see [I] for a pertinent review of the 
results along this line; also [TJ, [TT], [6]). Moreover, this kind of approach seems 
to be not able to account for partially opened cracks. In fact, in this case the 
transplacement field is continuous across the closed part of the crack although 
the material bonds are broken in the actual place. The description of partially 
closed cracks can have even stringent interest in the time-discretized proce- 
dure sketched above. In fact, when such a procedure is applied by taking into 
account a loading program described by time-dependent boundary conditions, 
a program implying even the growth of pre-existing cracks besides the nucle- 
ation of a new fracture at the n— th instant, it can happen that the cracks may 
close even partially, and then they re-open at subsequent time steps. Moreover, 
the discovery of physically appropriate crack patterns is another key point. 

Such questions have been tackled from a different point of view in [14] (see 
also additional comments in [15], [13]) by using tools from geometric analysis. 
The skeletal features of the theory are listed in the ensuing items. 

(z) Distinction is made between cracks and jump set of the transplacement 
field, as in [ID]. The latter set is however constrained to be contained within 
the crack pattern. Differently from all previous proposals, the crack pattern 
is described through measures over a fiber bundle with typical fiber the 
Grassmanian of 'planes' through B, the so-called curvature varifolds. 

(ii) A generalized notion of curvature can be associated with curvature var- 
ifolds. It enters the constitutive structure of the surface energy along the 
crack margins. The resulting energy differs from Griffith's proposal for the 
presence of the generalized curvature of the varifold and the energy along the 
tip in three dimensions. In this sense the model is an evolution of Griffith's 
scheme. 

The curvature-dependence of the surface energy has analytical advantages and 
permits the control of minimizing sequences. A theorem showing the existence 
of pairs of crack patterns and transplacement fields in appropriate measure 
and function spaces is then available (see [14J). It has some implications: 

(a) The crack pattern results a rectifiable set. Although it can be very irreg- 

of the material elements is then matter of modelling and can be even minimalist: 
in fact, one can choose to assign to every material element only the place that it 
occupies in the ambient space. I use to call Cauchy's bodies those bodies for which 
the minimalist approach summarized above is sufficient to represent the essential 
peculiarities of their morphology, and the representation of inner actions is just in 
terms of standard stresses. 
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ular, it has the features that our intuition assigns to a fracture. 
(6) Balance equations can be derived in weak form from the first variation 
even for a generic rectifiable set. Notice that if a crack is assumed to be 
coincident with the discontinuity set of the transplacement field only, to 
obtain balance equations, stronger regularity assumptions on the geometry 
of the crack pattern are necessary (see detailed analyses in [I], [5], [7], [TT] . 

ED- 

These consequences imply that, beyond the analytical advantages, the phys- 
ical meaning of items (i) and (ii) deserves additional analyses. They can be 
developed by going along the essential steps of the theory. 



2 

Consider the place B of a body, selected as a reference, as an open, connected 
set in the three-dimensional ambient space, with surface-like boundary ori- 
ented by the normal at each point, to within a finite number of corners and 
edges. If a 2D crack is formed in the deformed configuration u (B) - u the 
transplacement - and crosses a generic point u (x) with x in B, its 'direction' 
is locally described by tangent plane to the crack at u(x), when the crack 
is smooth. When the crack margins have a corner at u(x), a cone of planes 
has to be considered. Crack patterns can be however very irregular. One can 
accept a set as a representative of a crack pattern when it is just rectifiable as 
mentioned hitherto. So, an approximate notion of tangent plane is available 
in geometric measure theory (see [9]). 

Crack patterns can have a fictitious representation in the reference place B (let 
say through sets with zero volume measure) - the reference place is now mutant 
because the macroscopic structure of the body is changing with the nucleation 
and growth of a crack pattern. Let II indicate a two-dimensional plane or a 
straight line in three dimensional ambient space where B is contained. The 
pair (x, II) gives in B local information on the geometry of the crack crossing 
possibly u(x). Of course, up to when a real crack is not realized, any IT in 
the star at x can be a candidate to describe locally the direction of a possible 
crack pattern. The pair (x, II) can be viewed as a typical point of a fiber 
bundle Qk (B), k = 1,2, with natural projector n : Q\. (£>) — > B and typical 
fiber 7T~ 1 (x) = Qk,3, the Grassmanian of 2D— planes or straight lines associated 
with B. A k— varifold over B is a non-negative Radon measure V over the 
bundle Qk{B) (see [TJ, [2], [3], [TS], [IS])- The measure V has a projection 
over B obtained by using the projector of measures - indicated here by tt# 
- associated with the natural projection tv defining the fiber bundle Qk{B). 
Such a projection, namely n#V , is Radon measure over B and is also indicated 
for short by n v . It is called the weighed measure of the varifold and defines 
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the so-called mass M (V) of the varifold itself through the relation M (V) : = 
V (Gk(B)) — £v(^)- For the purpose of describing crack patterns through 
measures, the essential point is the possibility of constructing varifolds over a 
subset b of B. The subset b can have variegate nature. Here, for the physical 
purpose at hand, it is assumed that it is an admissible crack in the sense 
sketched above. Let 1-L k be the dimensional Hausdorff measure in M 3 , k = 
1,2. It is assumed that b is a H k — measurable, k— rectifiable subset of B. For 9 
a function in L 1 (bjH^, the approximate tangent k— spaed" 3 ! (here ID or 2D) 
T x b to b at almost every x in b is defined for 9H k \-b a.e. x G B. The symbol 
IT (x) indicates the orthogonal projection onto T x b. A varifold V bt g, restricted 
to b, can be then defined. It is called the rectifiable varifold associated with 
b, with density 9, and is such that 



Jo 



Gk{B) 



cp (x, n) dV b e (x,Tl)= ( 9 (x) ip (x, n) <m k 

Jb 



for any tp G C° (Qk(B)). Rectifiable sets can be considered a sort of generalized 
surfaces (see [I]). A subclass of them admits a notion of generalized mean 
curvature vector (see [2J, [3]). For members of such subclass (not all) a notion 
of second fundamental form can be defined (see [IB])- Here the attention is 
on varifolds admitting density 9 with integer values, the so-called integer 
rectifiable varifolds. They allow the definition of a special class of varifolds 
(see [19] ) which is essential for the ensuing development O- 

A new ingredient has to be inserted. It is a third-rank tensor field (x, IT) i — > 
A (x, n) G R 3 *®M 3 (g)R 3 * over Qk(B), with components Af. It plays the role of a 
generalized curvature. A varifold V is called a curvature A;— varifold with 
boundary if (i) V is an integer, rectifiable k— varifold V^e associated with the 
triple (b,9,H k ), (ii) there exists a function A G L 1 (Gk(B), M 3 * ® M 3 O M 3 *), 
and a vector Radon measure 5V such that, for every c/? G (^fc(jB)), one 
gets 

f (WD^ip + A l jD ul ip + Ajcp) dV (x, n) = - / dSV' (x, n) . 
Jg fc (B) v J J j J ' Jg k (B) 

The vector measure dV is called the varifold boundary measure [19]. The 
subclass of varifolds with generalized curvature A in L p {Qk{B)) is indicated 
here by CV k p (B). It is possible to show (see [18]) that if V = V hfi G CVj! (B), 
with p > k, V is locally the graph of a multivalued function of class C 1,a , 
a = 1 — ?, far from <9V^. 

• Varifolds with boundary can be used as descriptors of crack patterns: (?) The 
set b has the minimal geometrical properties of an admissible crack, at least 
in the sense mentioned above, (ii) The density 9 furnishes information on 



The choice k = 1, 2 allows one to treat in a unitary way surface and linear cracks. 
See [IB] , last chapter, for a nimble presentation of varifolds. 
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its possible faceted shapqjj. (in) The local orientation of the crack pattern 
is accounted for through II. (iv) Curvature is considered, although in the 
generalized (weak) form specified above, (v) The boundary of the crack - 
it includes the tip - is described by the boundary of the varifold. 

• Consider a smooth crack crossing a body in B and intersecting somewhere 
its boundary but maintaining the tip in the interior of B. A two-dimensional 
(k = 2) varifold V2 describes the crack, its boundary measure dVi is sup- 
ported by the entire boundary of the crack itself. To represent separately 
the crack tip, that is the part of the boundary of the crack in the interior of 
B, a specific one-dimensional varifold V\ has to be inserted. It is supported 
by the tip alone. Its boundary is supported by possible corners along the 
tip and the points determining the intersection of the tip with the external 
boundary of the body dB. The insertion of V\ allows one to assign later en- 
ergy to the tip of the crack. Different properties can be also assigned to the 
corners of the tip by using dV\. Of course, to capture the intuitive structure 
of the geometry under scrutiny, the varifolds V2 and V\ have to satisfy a cer- 
tain link. A definition presented in [H] specifies the link: a family {14}^ 
of A;— varifolds with boundary in d— dimensional ambient space is said to be 
stratified when tt# \dVk\ < f^Vu-i ^ or an ^' s - m ^ ne s P ec ial case treated here 
the condition of stratification reduces to 7r # \dV2\ < f/,y . 

• The choice k = 1,2 for constructing Q k (B) and the associated varifolds al- 
lows one to consider not only two-dimensional cracks with the relative tips 
but also additional linear defects (k = 1) which can be cracks included in 
very thin tubes - material bonds are broken along a line for some reason - or 
even dislocations. In the latter case, through a one-dimensional varifold one 
can describe dislocations emanating from a crack tip in a three-dimensional 
body. In the former case one manages crack patterns stratified over var- 
ious dimensions. Of course, the analogous description can be adopted in 
space dimension d > 2 and stratification of defects of various nature can be 
accounted for. 

A new form of the energy for a body undergoing fractures, based on the 
description of cracks in terms of varifolds, has been proposed in [T3] (see also 
comments in [15], [13]). Such an energy is indicated below by £ (u, {V*.} ,B). 
It differs from the expression (E (C, u) coming from the traditional Griffith's 
proposal and is an extension of it. For a three-dimensional body it reads: 



8 (u, {V k } , B) : = I e(x,u(x) ,Du(x)) dx + J2a k \ 

Jb ^ Jg k (B) 

+ j2(3 k M(V k )+ 1 M(dV 1 ), 

k=l 

5 If in a neighborhood of x there is a smooth surface, 6 = 1, when there is a net 
fold, 6 = 2, and so on. 
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where a>k, /3 k , 7 and pk are constitutive coefficients. In particular, a>k, (3 k , 
7 are positive numbers, so the contribution of the generalized curvature of 
the varifolds is always present, even if it can be extremely small. The density 
e (x, u, Du) is defined as the difference e (x, u, Du) : = e (x, Du)—w (u) between 
the bulk elastic energy e (x, Du) and the potential w (u) of external 
body forces. A number of comments on its physical meaning are necessary. 

• The addendum /3 2 M (V2) has the role of the last integral in (£(C,-u), that 
is the role of Griffith's surface energy: j3 2 has the same meaning of <fi in 
<E(C,u). Of course, the first addendum is the bulk elastic energy with density 
e (x, Du) - the body is then simple but anisotropic - and coincides with the 
first integral in <£(C,u). The other terms are not standard. 

• The addendum /3]M (Vi) counts energy along the tip. Such an energy is 
proportional to the length of the tip itself, namely to the mass M(Vi) of 
the one-dimensional varifold supported by the tip itself. 

• The term j~M.(dVi) adds possible energy concentrated at the tip corners 
where material bonds can be entangled in principle in a way different from 
the other parts of the tip. 

• The two addenda 



mark in a more pronounced manner the difference with respect to <E(C,u). 
They have pure configurational nature: they do not involve directly the gra- 
dient of deformation Du(x). The first one accounts for the (generalized) 
curvature of the varifold describing the surface of the crack, the second 
one includes the curvature of the varifold describing the tip. Influence of 
the curvature of the crack, above all in the proximity of the tip, has been 
recognized in [21] , a work devoted however to other aspects of fracture pro- 
cesses, namely Grinfeld's instability. The curvature here aims to account 
at macroscopic level of local microstructural effects at low scale, occurring 
in the cracking process. Curvature energy can be associated with bending 
effects in breaking material bonds when a crack is determined - consider for 
example a material in which the inner bonds are modeled through beam- 
like interactions. A point has to be stressed: Bending occurs in the current 
configuration u (B) while, as a function of x, A is defined over the ref- 
erence configuration B. However, when bending of material bonds in the 
actual configuration breaks the bonds themselves, such a bending has a 
configurational effect because it contributes to the mutation of B due to the 
nucleation and possible propagation of a crack. All configurational effects 
are measured in B, as it is commonly accepted. Additionally, analogies with 
different approaches to crack analyses in complex bodies (see [21]) and bod- 
ies with strain-gradient effects (see [25], [20]) can be called upon to enforce 
the interpretation of the presence of A as a configurational indicator of the 
effects due to latent microstructures. The terms including A in the energy 
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Fig. 1. Clamped two-dimensional sheet with a horizontal crack in the middle, sub- 
jected to Dirichlet boundary conditions - the arrows represent applied transplace- 
ments, the black zone is the clamping device. 

tell us essentially that, once boundary conditions are prescribed, one has 
to pay in energy for curving the crack. Consider a rectangular planar sheet 
of a homogeneous material half of it including a straight crack in the mid- 
dle as in Figure 1. Apply boundary conditions in terms of transplacement 
(Dirichlet boundary conditions) in mode 1. If the boundary conditions are 
such that the crack can growth, it remains straight, unless some additional 
agency occurs paying energy in curving it. 

• Consider a body in an initial configuration in which no crack is present. 
For example, B is a three-dimensional ball. After some loading program, a 
certain configuration in which a crack occurs is reached. Clearly one could 
consider such a configuration as the initial one. The possibility stresses 
the point that 'being cracked' is only a relative concept. Configurations - 
or better states - have to be compared to affirm that a body is cracked. 
Moreover, in sequences of configurations (or states) in which one excludes 
the possibility of restoring cracks by gluing the matter across crack fades, an 
order relation has to be considered. Such an order is given by monotonicity 
in crack patterns: if in a certain configuration there is a crack with respect to 
an uncracked configuration in the sequence, a subsequent configuration can 
have a crack pattern that coincides with or includes the previous crack. In 
terms of varifolds such a point of view is expressed by affirming the existence 
a family of comparison varifolds {Vfc} such that the family of varifolds {14} 
describing the actual crack pattern is constrained by \iy h < \i Vk for any 

k and V k G CV£ k (£>). The assignment of ji4j does not mean that one 
is considering in a given configuration a pre-existing crack pattern always, 
because the comparison varifold family can be also empty. 

• The assumption used so far that the ambient space is three-dimensional has 
been accepted only to be close to the standard physical intuition. There is no 
obstruction to consider higher dimensional spaces. Consider the dimension 
of the ambient space to be d > 2. The previous treatment can be extended 
to this case straight away. The only modification, including the expression 
of the energy is that the range of k must be considered tobe 1 < k < d — 1, 
at least in principle. The rest remains unchanged. This choice allows one to 
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consider stratified families of defects at various scales with various physical 
meanings, depending on the circumstance. Of course, there could be cases 
in which k could not take all the natural values from 1 to d — 1. Even in 
three-dimensional ambient space one could have only k — 1, for example, 
being in the condition to consider only linear defects like nets of disloca- 
tions. Such would deserve perhaps a treatment a part. The case of 
four manifolds, d = 4, could be called upon when the description of crack 
in relativistic elastic bodies is considered, with all the necessary changes in 
the representation of the energy, adapted to the relativistic setting. In this 
case, however, care must be taken because in the relativistic setting a rep- 
resentation of continuous bodies based on the back-to-label representation 
seems to be preferable. The circumstance then would change the stage and 
would include a sort of 'mixed' representation of the energy. 

In the setting described so far, a minimality requirement is prescribed for the 
energy: 



Minimize E (u, {V^} , B) with V k in CV£ k (£>), comparison varifolds jV^j, u 

in an appropriate function space, with assigned boundary conditions. 

Solution to this problem, if any, is a pair (u, {14}). Minimization over a class 
of varifolds has the meaning of minimization over a class of bodies: every pos- 
sible crack pattern represented over B (remind that nucleation and\or growth 
of cracks occur in the current place) defines a body. Different crack patterns 
indicate different bodies: B, in fact, changes. The family of varifolds rep- 
resents the crack pattern - so it selects a body - and the field x i — > u (x) 
describes the deformation of such a body. The two 'objects' are correlated. In 
fact, in the actual configuration possible nucleation, growth and/or opening of 
a crack are consequences of the deformation, the varifolds over B are represen- 
tatives of what happens in the actual configuration. The choice of the function 
space hosting the generic u is then another key point of the treatment. It has 
to be linked to the varifolds supported on B. 



In continuum mechanics the choice of function spaces as ambient for solutions 
to equilibrium or evolution problems has constitutive nature. The properties 
of the members of a given space carry a physical meaning about. The charac- 
teristic physical features of the problem under analysis have then to address 
the functional choice. In the case under scrutiny, the idea is that the material 
is elastic-brittle. It means that the material is elastic up to a certain threshold 
after which a crack is created while outside the crack the material is still in 
elastic phase. The threshold can be expressed in terms of deformation, stress, 
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energy, depending on circumstances. Moreover, as it will be clear later, the 
threshold can be also not expressed directly. It is, in a sense, included in the 
choice of the function space that one selects. 

Let us focus the attention on pure elasticity first. If the material is purely 
elastic, it can deform at will, without any threshold and, in principle, with- 
out end. The deformation can be also perfectly recovered, after unloading. In 
this sense, phenomena like cavitation in solids are ascribed to elastic-brittle 
behavior rather than perfect elastic setting. If this view is accepted, the conse- 
quence would be that a perfectly (hyper)elastic body should be such that any 
compatible transplacement (or displacement, depending on the choice) field 
does not describe the nucleation of fractures or holes. 

It is well known that under conditions of polyconvexity of the elastic energy 
density, minimizers of the elastic energy of a simple body can be found in the 
Sobolev space W 1,p {B, IR 3 ) - a space hosting maps with first distributional 
derivative having integrable p— power, i.e. the first derivative is in L p . However, 
when p < 3 non negligible is the presence of transplacement maps with graphs 
admitting boundaries with projections into the interior of B. Such boundaries 
describe the formation of 'holes' and/or open 'fractures' of various nature, so 
they are undesirable when a purely elastic material is under analysis, at least 
if the view on elasticity sketched above is accepted. 

The difficulty can be overcame. In fact, there is a global way to check - even- 
tually to control - the presence of boundaries in the graph of a map with pro- 
jection into the domain of the map itself through linear functionals. Indicate 
them by G u , with the indices u suggesting that the functional G is associated 
with the transplacement u. In the case treated here, such functionals are linear 
over smooth 3— forms with compact support inBxl 3 . The general technique, 
which is valid in any finite dimension d of the ambient space, is described in 
the monograph [17]. Here, I sketch only minimal ideas, furnishing the essential 
picture and physical interpretations in the case under scrutiny. 

Preliminarily, remind that a r— vector over a linear space E is a rank— r skew- 



symmetric tensor, that is an element of the skew-symmetrization of E®...®E. 
The space of r— vectors is indicated here by A r (E). It has a natural dual 
A r (E). Any map of the type u : B — > A r (E) is called a r— form0- The space 
of all r— forms of the type just defined is indicated by V r (E) 

For the mechanics treated here interesting is the case of 3— vectors over IR 3 x 
IR 3 , the space hosting the graph of the deformation. Consider a deformation 



6 More specifically, one should consider £ as a real vector bundle over B of fiber 
dimension d, that is a family of d— dimensional vector spaces parametrized by points 



r 




of B. 
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(transplacement) x i — > u (x) G M 3 , x G i3. Amid all possible 3— vectors over 
IR 3 x M 3 , it is possible to define at every x G B the 3— vector M (Du) asso- 
ciated with the gradient of deformation Du. Detailed definition and various 
properties can be found in [17] . here, for the expository purposes declared at 
the beginning, it is only necessary to know that at each x its components 
are the entries of Du (x), adjDu (x), detDu (x). In M (Du), then, all elements 
characterizing the deformation of lines, areas, and volume of the body in B 
are included. In this sense, M (Du) characterizes completely the deformation. 
Its dual counterpart - the value at the same x of some form in T> 3 (K 3 x M 3 ) 
- is then a sort of generalized stress. 

Given a transplacement x i — > u (x) G M 3 , x G B, the 3— current integration 
G u (current for short) over the graph of u is defined to be a linear functional 
over smooth 3— forms u G T> 3 (M 3 x IR 3 ) with compact support in B x R 3 , 
namely 

G u (u) : = / (uj (x,u(x)) , M (Du(x))) dx, 
Jb 

where the angle brackets indicate the natural action over M (Du) of its dual 
counterpart] 7 [ Essentially, G u (u) plays the role of generalized internal power. 
The number M (G u ) indicates here the so-called mass of the current and 
is defined by 

M(G U ) := / \M(Du(x))\ dx, 
Jb 

where \M (Du(x))\ is the modulus of M(Du(x)), evaluated in the standard 
way for tensors. The symbol \G U \ indicates the total variation of the current 
and is defined as usual for functionals. A boundary current can be associated 
with G u : it is indicated by dG u and defined by duality, that id" 8 ! 

dG u (w) ■= G u (du) , Vw G V 2 (B x M 3 ), 

with T> 2 (B x M 3 ) the space of 2— forms over M 3 x M 3 with compact support in 
B x M. 3 . The notion of boundary current has not only formal nature. It has an 
immediate physical interpretation: when the graph of u is free of boundaries 
inside the interior of B, dG u (oo) = for any u G T> 2 (B x M. 3 ). Essentially, 
this zero boundary condition prevents the formation of cracks or holes inside 
the actual place u (B) of the body. Such a condition has been used (see [T7] 
and the other references of its authors mentioned therein) to define a class of 
transplacements - the so-called weak diffeomorphisms - which is 'consti- 



7 A bit more precisely, the current is defined by taking the rectifiable part of the 
graph of u, that is the part of the graph that can be seen as the graph of Lipshitz 
maps. 

8 The so-called external differentiation over forms is indicated by d and acts as 
d:V n (E) ->V n+1 (E). 

9 The definitions of currents and related boundaries can be also available in spaces 
with higher dimension (see [UJ for the complete theory). 
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tutively' an appropriate choice for describing what one imagines to be a pure 
elastic deformation, as sketched above. 

Here the situation is a bit more complicated. The deformation has to be a weak 
diffeomorphism outside a subset of the support of the varifolds describing the 
crack pattern, if the minimizing procedure provides a non-empty family of 
minimizing varifolds. Inside that subset, the transplacement admits jumps, so 
it is not purely a weak diffeomorphism. 

Extended weak diffeomorphisms are then necessary. They must have the phys- 
ical properties just indicated. Such properties can be summarized in a formal 
definition. 

Definition 1 Assigned a stratified curvature varifold V = {Vfc}^~| with bound- 
ary, i.e., Vk G CV Pk , a map x i — > u is said to be an extended weak diffeo- 
morphism (in short u G dif 1 ' 1 (B, V, M 3 ) ), when 

(i) mGL 1 (B) and is a.e. approximately differentiable, 

(ii) \M{Du)\ G L 1 (B), 

(Hi) det Du (x) > for almost every x G B, 
(iv) for any f G C c °°(i3 x R 3 ) 



/ / (x,u (a?)) det Du (x) dx < sup fix, w) dw, 

2 

(v) 7r# \dG u \ < Hy u + 7r# \dV\\ as measures on B. 
i=i 

The definition has natural extension in ~R d : the summation in the item (v) 
should be extended up to d—1. Of course, the definition of currents and related 
boundaries holds in dimension d: in that case M (Du) is tested over d— forms. 
Here and in the whole paper, the restriction d = 3 is essentially motivated 
by the physics under scrutiny (see [H] for the abstract theory). The first 
item indicates the possibility of evaluating the gradient of deformation. The 
second item is another regularity condition. It implies that one can in principle 
measure the average of the gradient of deformation, the volume change, the 
overall deformation of surfaces. The third item is the standard condition that 
a transplacement be an orientation preserving map. Item (iv) is the condition 
mentioned at the beginning of Section 1. It permits to move B along u into 
a region u (B) in such a way that self-contact between parts of the boundary 
dB be allowed while self-penetration excluded. Notice that in item (v) the 
action of the projector 7r# on the total variation of the boundary current 
is motivated by the fact that the latter behaves substantially as a measure. 
Essential properties for the space of extended weak diffeomorphisms are shown 
in [T5] (see there the relevant theorems and proofs). 
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Standard weak diffeomorphisms are W 1,1 [B, IR 3 ) maps that satisfy the items 
(ii), (Hi), (iv) in previous definition while item (v) which is substituted by 
the zero boundary condition dG u (u>) = for any u G V> 2 (B x ]R 3 ). From a 
kinematic point of view, in going from dif 1,1 (B, IR 3 ) to dif 1 ' l (B J V,M?) ) one 
transits from pure elastic setting to elastic-brittle behavior. 

In this sense, the requirement of minimality of the energy includes the possi- 
bility of finding minimizers - if any - in terms of weak diffeomorphisms and 
null varifolds, and in terms of extended weak diffeomorphisms and non-null 
varifolds. The transition from a situation to another is morally the threshold 
from the elastic to the elastic-brittle behavior. In this sense, also, there is 
no need in principle of adding another condition defining the threshold itself. 
There could be also minimizers for which the transplacement field is simply a 
weak diffeomorphism but the varifolds are not null. This situation describes 
presence of closed cracks only: in other words, material bonds are broken but 
the crack remains closed and the transplacement field does not jump across 
the crack facies. Such a situation can occur in a step-by-step minimization 
program obtained by updating in time steps the boundary conditions and 
requiring minimality of the energy at each step. At the step n — 1 a crack 
pattern can occur, at the step n the deformation closes the cracks, at the step 
n + 1 there is a purely elastic continuation, then, at further steps, new crack 
patterns accrue. 

Of course, proving existence of minimizers is a crucial step for attributing sense 
to the previous reasonings. Existence depends on the characteristic properties 
of the energy and the boundary conditions. Once the existence of minimiz- 
ers is established, the characterization of them along the physical suggestions 
collected above is matter of regularity theorems. The question is open. Actu- 
ally, any regularity theorem is available. Different is the case of the existence 
problem. 
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The existence result for the minimum problem stated above for the energy 
£ (u, (Vfc},£>) of an elastic brittle solids has been proven in [TJ]. Boundary 
conditions of Dirichlet type can be presumed. They are given by prescribing 
the transplacement field along the boundary dB of B. 

The discussion of the existence is developed by taking first a subspace of 
dif 1 ' 1 ^, V,R 3 ), precisely the space dif p '\B,V,M. 3 ) defined by 

di f p,1 {B, V, M 3 ) := jw G di f 1,l [B, V, IR 3 ) | \M (Du)\ G LP (£>)} , 
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for some p > 1. Essentially, the choice of dif p ' 1 (B, V, M 3 ) with p > 1 is a re- 
quest of additional regularity which is sometimes necessary for physical needs. 
Combination with the space of varifolds allows one to recognize a natural am- 
bient in which the existence of minimizers of the energy £ (u, {V^} , B) can be 
investigated. Such a space is indicated by •A qpK |y i (£>) and defined by 

\ P ,K,{v k } ( B ) : = 04)) I ^ec^^^erty,! 3 ), 

{Vfc} is stratified, \\u\\ Loo{B) < K,fi Vk < v Vk ,Vk = 1,2} , 

where V\ and V2 are comparison varifolds describing possible initial cracks. In 
particular, the subspace 

A Z,K,{V k ) ^ ■= W) G \p,K,{^} (B) 1 U (X) = U ^ X ^ Xe dB »> } 

with dB u the part of the boundary of the body where the transplacement field 
is prescribed, takes into account the boundary conditions of Dirichlet type 
mentioned above. 

In all these definitions, another regularity requirement is prescribed. In fact, 
the condition IMI^o^g) < K imposes that the essential supremum of u is 
almost everywhere - with respect to the Lebesgue measure - bounded. In 
fact, a priori it is not possible to exclude that, if one is able to prove under 
some conditions the minimality of the energy E (u, , B) over some space 
of extended weak diffeomorphisms and varifolds, the minimizing varifold does 
not describe a fragmentation of the body - let say a crack cutting a piece 
of matter from the rest. In this transplacement field could be such 

that the cut piece, now free from boundary conditions, can be translated 
rigidly to infinity. By imposing that IHI^oorg) < K, with K a real number, 
then, one wants to avoid the situation just sketched. So, in this sense the 
assignment of K has not properly constitutive nature. It is not related to 
some property of the material, rather it is a parameter selecting admissible 
deformation processes, admissibility considered with reference to the possible 
unconstrained extraction of pieces of matter from the body. 

For d the dimension of the ambient space and k ranging from 1 to d — 1, 
analogous definitions of A u ° r , r ,-, 1 (B) hold and the theory can be generalized 

& q,p,K,\V k j v > J a 

(see relevant results in [2]). 

Another crucial point in the path leading to the proof of existence theorem of 
minimizers of the energy is the discussion of the structural properties of the 
energy. They have constitutive nature, of course. 

In non-linear elasticity of simple Cauchy's bodies, common assumptions about 
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the structure of the energy e (x, u (x) , Du (x)) are well known. By indicating 
by M^ 3 the space of 3 x 3 matrices with positive determinant, the energy 
density e is considered as a map 

e : B x R 3 x M 3 + x3 ->- [0, +oo] 

with values e (x, u (x) , Du (x)). Remind that positiveness of the determinant 
of Du (x) is the condition assuring that the transplacement be orientation 
preserving. The properties H1-H4 below are then assumed to hold. 

HI e : B x R 3 x M 3 + x3 — > [0, +oo] is continuous in (x,u). 
H2 The map Du(x) i — > e (x, u (x) , Du (x)) is polyconvex: that is there 
exists a Borel function Pe acting as 

Pe : B x M 3 x A 3 (M 3 x ]R 3 ) ->• R+ 

with values Pe (x,u(x) , £ {x)), which is continuous in (x, u) for every £ G 
A3(M 3 x M 3 ), convex and lower semicontinuous in £ for every (x, m), and such 



that! 10 1 Pe (x, u, M {Du)) = e (x, u, Du) for any list of entries (x, u, Du) G 
B x R 3 x M 3 + x3 with det Du > 0. 
H3 The energy density e satisfies the growth condition 

e(x,u,Du) > Ci |M (Du)| r . 

H4 For every x G S and D« G M 3 + x3 , if for some u G M 3 the inequality 
e (x, u, Dm) < +oo is satisfied, then det Du > 0. 

The physical nature of these assumptions is discussed in various treatises (see 
|22j . |23j). The standard presence in the polyconvex energy of the determinant 
of the gradient of deformation and the relevant adjugate is summarized here 
in the functional dependence on M (Du). The choice is not only formal. It 
furnishes a rapid path toward the extension of the treatment to d— dimensional 
cases (see [T7]). 

However, the essential point is to underline that, in the setting explored here, 
H1-H4 do not need to be supplemented by additional structural assumptions 
on the energy to assure the existence of minimizers of £ (u, {Vk} , B). The 
relevant theorem reads as follows: 

Theorem 2 (|14|) Assume K > 0, q,pk > 1, and T4 G CV£ k (B) for any k. 

If there exists (u , {Vg}) G A u ° pR r-, (B) such that £ (u , {Vg} , B) < +oo, 

then £ (u, {Vk} ,B) attains in that space the minimum value. 
Proof is presented in [H]. Here just comments have to be added. 
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The dependence of u and Du on x is now suppressed for the sake of brevity. 
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Since no additional structural hypotheses besides H1-H4 of standard non- 
linear elasticity need to be added, the information about the possible nu- 
cleation of a crack or growth of an existing one is furnished by the presence 
of the terms ruled by varifolds in the energy and the constitutive choice of 
A u ° „ r - (B) as functional setting. It is just the latter choice that avoids 

q,p,K,{V k \ v ' 

the introduction of an external criterion for the nucleation of a crack or the 
growth of an existing one. In fact, the energy is minimized over a class of 
possible bodies. 

More in general than other descriptions, it is possible to determine the weak 
form of balance equations for crack patterns which are just rectifiable sets. 
The result is not discussed here for the sake of conciseness. It is presented in 
[H] and opens the way to computational opportunities not explored yet. The 
balance equations derived naturally in [13] for very general crack geometries 
- as mentioned above the crack pattern has to be just a rectifiable set - 
are the ones obtained by horizontal variations, that are variations of the 
reference placd~^~l. Thus they have configurational nature: they involve in 
fact the Hamilton-Eshelby tensor and non-standard terms deriving from 
the variations of the terms including the varifolds, the ones directly related 
with the geometry of the crack pattern. 

The existence result holds also at dimensions greater than 3, provided that 
the obvious variations in previous definitions (see [14] , [To] for the abstract 
theory) . 

An analogous result holds also for the generalized energy 

„ d— 1 „ 

£ (u, {V k } ,B):= / e(x,u(x) , Du (x)) dx + V a k \ (p k ( Aha ) dV k + 
+ J2f3 k M(V k )+jM (dVt), 

k=l 

with <j) k (•) a convex real-valued function satisfying the condition <p k (t) > 
ct Pk . The interest of this remark is not only technical. It gives a grater degree 
of freedom in selecting further constitutive structures under the suggestions 
of possible experimental evidences and numerical tests. 
The existence result does not exclude the possibility that the minimization 
procedure foresees stratified varifolds supported on the boundary of B. In 
this case one can say that a boundary crack appears. The meaning of such a 



11 See [8] for clear explanations on the connection between horizontal variations de- 
scribing the potential movement of defects and the balance of configurational forces. 
Take into account also that there is basic difference between the balance of config- 
urational actions associated with macroscopic mutations of the reference place and 
the balances of standard actions generated by the deformation. In the conservative 
case and with reference to smooth fields, the former balances are essentially the 
pull-back in the reference place of the latter balances. In general it is not so. The 
difference has been evidenced first in [12] vol.1, pages 152-153. 
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boundary cracks is not exotic. In fact, on a part of the boundary where the 
transplacement is imposed a crack can occur so that the boundary condition 
is 'broken'. Consider for example a beam jointed at one of its ends. Boundary 
conditions and properties of the material can be such that a crack occurs 
just at the interface between the joint and the beam. However, in principle a 
boundary crack can appear on a free part of the boundary. Such a situation 
can describe the fragmentation of a thin film at the boundary, which is, 
essentially, the abrasion of the boundary itself. 
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The technique discussed previously is a general tool for the description of phe- 
nomena in which energy can in principle be concentrated over submanifolds 
of a certain manifold, and this energy depends on the geometry of the sub- 
manifold itself. It can be used to analyze either specific situations of physical 
interest or to formulate and analyze abstract mathematical problems. 

• Phenomena of physical interest that can be described by using the tools 
mentioned hitherto deal for example with the mechanics of linear defects 
like dislocations and discontinuity surfaces. For dislocations the choice of 
the varifolds play a crucial role. For interfaces, the functional setting has 
to be changed and a new special class of extended weak diffeomorphisms 
arises. The new choice requires proof of completeness of the new functional 
class and evaluation of the applicability of lower semicontinuity results. 

• Another point of discussion is also the evaluation of the crack nucleation 
and growth in complex bodies. The adjective 'complex' distinguishes bod- 
ies characterized by a prominent influence of changes in material texture 
(the microstructure) on the macroscopic behavior, an influence exerted 
through inner actions requiring a representation going beyond the common 
picture in terms of standard stresses. Quasicrystals, ferroelectrics, magne- 
toelastic materials, polymeric bodies of various nature, including elastomers, 
fullerene-based composites, porous bodies, bodies with continuous distri- 
butions of dislocations, multiphase materials are paradigmatic examples. 
Notwithstanding the variety of special models, a unitary picture of the me- 
chanics of complex bodies exists. Within it, the representation of bodies 
goes beyond standard Cauchy's approach in the sense that every material 
element is viewed as a system rather than a black box individuated by a sin- 
gle point in the ambient space, which is Cauchy's view. Such a description is 
multifield and intrinsically multiscale. A morphological descriptor field 
x i — > v (x), x e B, of the essential geometrical features of the material mi- 
crostructure is then introduced. To construct the essential structures of the 
relevant mechanics, it is just necessary to presume that v (x) is an element 
of a set M which has just the structure of a differentiable manifold. It is 
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assumed to be finite-dimensional for the sake of simplicity. M. is called the 
manifold of substructural shapes. The interest of mentioning here com- 
plex bodies is motivated by data showing that the microstructural changes 
may influence in non negligible way the force driving crack tips along evolu- 
tion processes. Theoretical analysis of this phenomenon within the setting 
of the general model building framework of the mechanics of complex bodies 
(that is without specifying the type of microstructure) has been developed 
in |2lJ~^l from a point of view different from the one adopted here. The 
simplest extension of the theory discussed here to complex bodies is given 
by an energy of the type 

£ (it, {Vk} . B) : = / e (x, u (x) , v (x) , Du (x) , Du (x)) dx + 
Jb 

2 „ 2 

dV* + £&MM+ 7 M (dVt), 

k=l k=l 

with the natural modifications in generic dimension d. In analyzing the ex- 
istence of minimizers, an essential point is the choice of the space hosting 
the morphological descriptor maps. Such a choice could require the embed- 
ding of M. into a linear space. Such embedding always exists because M. is 
finite dimensional, also it can be isometric when Ai is Riemannian. In all 
cases, however, it is not unique so that it becomes a ingredient of the model, 
a sort of constitutive choice. Another point is the link of the jump set of 
the morphological descriptor field with the varifolds. Here the underlying 
physics is subtle. In principle one can accept that x i — > v (x) may have 
jumps even outside the support of the varifolds and there it may be even 
continuous. Jumps outside the varifolds can be justified by the formation 
of domains of microstructures, like polarization or magnetization domains. 
The meaning of the possible continuity on the support of the varifolds is 
associated with the question whether in cracking a body one alters along 
the margins of the crack the microstructure, in a sense determining a new 
type of microstructure, or, else, the microstructure remains the same across 
the margins of the crack. In the philosophy of continuum mechanics, the re- 
mark above coincides with asking whether a crack just divides neighboring 
material elements or breaks the material elements met in front of the tip. 
The answer cannot be definitive and is matter of modelling. The situation 
becomes also more complicated when the structure of the surface energy 
involving the generalized curvature of the varifold is enriched by making 
more articulated assumptions. Relevant investigations are actually open. 
• A point which may deserve to be noted is that the scheme discussed in 
previous sections has intrinsic similarity with the general framework of the 
mechanics of complex bodies which has been sketched rapidly in the last 



12 There one can find appropriate references to works presenting and discussing 
experimental data. 
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item. In fact, when curvature k— varifolds with boundary are chosen to rep- 
resent cracks, in principle the region where they localize is not known - 
in other words one does not know where the support of the varifolds is 
placed in B, that is where the crack is. In this setting, every point can be 
crossed in principle by a crack. The minimization procedure tells us that the 
crack is here or there, before nothing is known about its position. Instead 
of assigning to each material element a morphological descriptor v selected 
in a finite dimensional differentiable manifold Ai, one is then assigning to 
each material element a measure. In this sense, the scheme discussed here is 
driven by the ideas of the mechanics of complex bodies, and, in some sense, 
it goes beyond them a bit. 
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